Externally-driven collisions of domain walls in bistable systems near criticality 
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Multi-domain solutions to the time-dependent Ginzburg-Landau equation in presence of an ex- 
ternal field are analyzed using the Hirota bilinearization method. Domain-wall collisions are studied 
in detail considering different regimes of the critical parameter. I show the dynamics of the Ising 
and Bloch domain walls of the Ginzburg-Landau equation in the bistable regime to be similar to 
that of the Landau-Lifshitz domain walls. Domain-wall reflections lead to the appearance of bubble 
and pattern structures. Above the Bloch-Ising transition point, spatial structures are determined 
by the collisions of fronts propagating into an unstable state. Mutual annihilation of such fronts is 
described. 
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I. INTRODUCTION 

Localized and labyrinth structures in excitable media 
are considered to be used as carriers of binary infor- 
mation for specific applications [H, [!]• Such structures 
appear in monostable and bistable regimes of magnetic, 
electrical, chemical, and optical systems. In particular, 
different realizations of memory and logic elements with 
bubble-forming bistable media seem to be technologically 
promising because of a high stability of the bubble struc- 
tures. Prototype bubble-forming systems are 2D ferro- 
magnets including ferrofluids 0, 0j , and 2D ferroelectrics 
including electrically active hquid crystals 0,0 . In criti- 
cal parameter ranges, these systems create pattern struc- 
tures (lamellae, labyrinths) which are not periodic [3, 0|- 
To be precise, I mention that experimental observation of 
bubbles and lamellar patterns in solid ferroelectric films 
does not manage (unlike in liquid crystals where electro- 
elastic coupling is weaker) , although a 2D bubble like 
structure has been seen in a system of a finite-size geom- 
etry Unlike patterns in oscillatory critical systems, 
the above structures are built of domain walls (DWs) con- 
necting stable states of opposite orientation of an order 
parameter. Similar patterns and bubbles are observed 
in excitable (reaction-diffusion) chemical systems below 
the critical (bifurcation) point [l^, [lH . It is hoped that 
difficulties of early magnetic and dipolar bubble tech- 
nologies reviewed in (l^nisj can, nowadays, be overcome 
thanks to the advance of miniaturization. One observes 
increased interest in studies of DW complexes with rel- 
evance to data storage and processing devices which are 
based on ferromagnetic and ferroelectric nanoelements. 
In particular, nanowire structures are under intense in- 
vestigation because of t he p ossibility of the most dense 
information packing [l^ [TBI . The importance of critical 
effects in nanosystems grows with decreasing in their di- 
ameters because of a significant decrease of the critical 
(Neel) temperature [T6l.[l7j. 
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Generic critical properties of the magnetic and electric 
media are described with the time-dependent Ginzburg- 
Landau equation which enables qualitative analysis of 
DW behaviors in all bistable systems. The DW solu- 
tions to the ID Ginzburg-Landau equation are of Ising 
(named Neel DWs also) or Bloch type. Under external 
stimulation the DWs move and their collision properties 
determine morphology of resulting pattern structures. In 
the present paper, external-field-induced DW collisions 
in ID Ginzburg-Landau systems are studied with refer- 
ence to bubble formation. I show similarities of the pro- 
cess to the bubble formation in ferromagnetic (or ferro- 
electric) wire far from the criticality (a system described 
with the Landau-Lifshitz-Gilbert equation), referring to 
my previous study of the problem [18] as I throughout 
the text jl^] . The creation of bubbles and DW patterns 
is described with connection to the property of elastic 
refiection of DWs observed in excitable media. How- 
ever, I emphasize that the present study is not related 
to widely investigated pulses in monostable regimes of 
excitable systems [l^ [lo] . 

In sections II and III, the field-induced collisions of 
DWs and phase-fronts, respectively, are studied. In sec- 
tion IV, the relevance of predictions on DW and front 
collisions to pattern formation is discussed. 



II. DOMAIN WALL COLLISIONS 

Let us consider ID Ginzburg-Landau equation for the 
ferromagnetic (ferroelectric) wire in a (longitudinal) ex- 
ternal field directed along the spontaneous magnetization 
(polarization) 

dm d'^m n * i i2 tt 

a-^ ^ J -^-^ + Pim + -fi\m\m + -fH. (1) 

Here m denotes a complex order parameter (a two- 
component magnetization or polarization), /3i determines 
the distance from the criticality {/3i = —(32 at the phase- 
transition point), H denotes an external (magnetic or 
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electric) field intensity (it can take positive or nega- 
tive real values). Basing on arguments related to the 
time-reversal symmetry, I show that the described with 
(HI) DWs have to reflect upon their collision in the pa- 
rameter region where they form magnetic bubbles. Un- 
der the assumption <^ (/3i -I- 132)^^'^ / fJ-^^'^ , the do- 
mains of the stable phase of the system correspond to 
m ~ ±-\//3i +132/ y/V^, and the single- DW solutions to (P) 
take the form 



/3i -f 1 - e^" , .^^ /?i-3/?2 e" 
m=\l z — + 1O2 



77 — k{x — Xq) — 



^ 1 e^") y ^ 1 -I- e^'' ' 



(2) 



where 



and 



2/32 



\k\ = \l^ for 9^±1, 



1.1 = ./^ for 



= 0. 



(3a) 



(3b) 



These equations describe the Bloch DW and the Ising 
DW, respectively. Only in the limiting case H = 0, the 
above solution is strict while, for _ff 7^ 0, it neglects a 
small field-induced deformation of the DW profile (an 
space asymmetry arising with respect to the DW cen- 
ter; details are given in Appendix A). The DWs move 
under the condition TJ 7^ in the direction dependent 
on the sign of k. It is seen from (l^. pap that the Bloch 
DWs can exist only in the regime /3i > 3(32; therefore, 
at /?! — 3/3^, the Bloch DWs transform into the Ising 
DWs which are the only known stationary front solu- 
tions in the regime —(32 < (3i < SB-? |2l|. Stationary 
two-DW solutions to ^ are given in In absence of 
the external field H — 0, the double- wall solutions de- 
scribe localized static bubbles built of pairs of different 
type (Ising and Bloch) DWs. The size of the bubble can 
be changed by the application of a field H ^ because 
of inducing the motion of its fronts (DWs) and, eventu- 
ally, their collisions. The collision result (DW reflection 
or annihilation) determines final spatial structure created 
by DWs, similar to whose of ID ferromagnets far from 
the critical regime (I). In the present work, I omit the 
discussion of the formation of bubbles built of two Bloch 
DWs or two Ising DWs mentioning previous studies of 
Refs [2^. Such objects are unstable for H — (initially 
stationary DWs repeal or attract each other), thus, they 
are less important for potential applications than Ising- 
Bloch bubbles (although, the stability of the Ising-Bloch 
bubbles against temporal external perturbations is a sep- 
arate question discussed in [24]). 

Obviously, Eq. ([T]) is not time-reversal invariant, while 
the reversibility of a dynamical system is a crucial prop- 
erty in terms of the stability of any localized structure 
(e.g. a bubble) T,'2^. Due to this irreversibility, solitary- 
wave solutions to ([T|) are, in the general case, relevant 



only to the limit of large positive values of time (t — ^ 00), 
while they are irrelevant to the opposite limit {t — ;> —00), 
in particular, because of infinite growth of the energy 
with t ^ —00 which is indicated by nonzero values of 
the dissipative function (I). With relevance to studying 
the DW collisions, and thus to investigating the stability 
of DW complexes, I analyze the dynamics of the bubble 
in the limit of large negative values of time (in this limit, 
the colliding DWs are noninteracting objects), modifying 
the evolution equation ([Ij by reversing the arrow of time. 
The inverse-evolution equation takes the form 



drfi 



J 



+ (3im + f32ih* — /i|m| fh — jH. (4) 



Let me mention an analogy to the necessity of doubling 
the number of degrees of freedom when formulating the 
dynamics of dissipative (classical or quantum) systems 
within formalisms relevant to the whole length of the time 
axis ^26j, J^, (see also I), in particular, the dynamics of 
essentially dissipative (reaction-diffusion) systems whose 
excitations are overdamped f28| . 

Applying the Hirota's bilinearization method, follow- 
ing the substitutions m = (71//1, m = — (?2//2 in (HI) 
and ([4]), where /i, /2 take real values, one arrives at the 
secondary equations of motion 



i-aDt + JDl)g, ■ f, + - X)g,h 

+p29lh+lHfl^0 

JDl)g2 ■ /2 + Wi - A)52/2 

+/32.g2*/2 + 7^/1 = 

JDlf2- f2+ (^929*2 -Xf^^O 



{aDt 



(5a) 



(5b) 



Here, 



D^-D^bix, t) ■ c{x, t) = [d/dt - d/dt )™ 

X {d/dx ~ d/dxYh{x, t)c{x , t)\,^,'^t^t' . 

The above breaking of (H]), (U)) into the pairs the sec- 
ondary equations ([5a| . (I5b[) . respectively, is nonunique 
but it leads to the equations of the lowest possible order 
in 5i(2)i /i(2) (bilinear ones). Upon the inversion of the 
arrow of time {t — >■ — t), the secondary dynamical vari- 
ables transform following -J> -52(1), /i(2) ^ /2(i) 
while ([5a|) transform into (I5bp and vice versa. 

For \ = (3i+ (32 and -fH < (/3i + Bifl'^jii^l'^, (the 
field is much weaker then the coercivity value; thus, m w 
±-\/ f3\ + (32I ^fji inside the domains and a deformation of 
the DWs induced by H is negligible) , following Appendix 
A, I find a two-DW (bubble) solutions in the form 



gi = a(l - j/e 



2»)1 



.2)): 



2»72 



g2r;i+2»;2 
2))i+2jj2 



l/2^')l 



(1 



:.2l)s 



/i = I -I- vef-'^ -I- i/e"''" -f e 
r\j — kj{x — xoj) — ujjt, 

92 = a{l - i^c^''! - ve^^^ + fflm+^fi^^^ _^ i6z/i/2gi)i _ g: 



(6a) 



2i)2 



,2»?i 



„2i)2 



/2 = I -h ve~ '^ + ve"'' + e 
fjj — kj{x — xoj) + ujjt, 



2fji+2fj2 



(6b) 
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where 



//3i-3/32 



fci = T 



2J 



/?i+fe-2V/?2(/3i+/32) 
/3i+/32 + 2V/32(/3i+/32)' 



(7) 



In the Hmit iJ — > 0, 51 = 52, /i = /2, and both the fields 
m and — m coincide with the static bubble solution to ([T]) 
written in [2^. 

Studying the DW collision, I analyze the limit f — > 00 
of the double DW solutions to ([T]), the formula (l6al) . and 
the t — > —00 limit of the solution to (|1]), the formula (j6bp . 
in the closest vicinity of the i-th DW center, noticing that 



]^ _ j,2j7,+ln(i/) gr,,+ln(i/)/2 
a- ^ — -r^ + i*i2^^:;— — ^„.^,„,,a (8) 



for : 



for 



77, + ln(i^)/2 > -qk + ln(i/)/2 and i ^ fc. 



1 



„2,ji+ln(iy) 



' I e2))i+ln(i/) 



i(5i2& 



gf)i+ln(i.)/2 
1 -|_ Q2f,i+\n(v) 



(9) 



% + ln(i/)/2 > ryfc + ln(i')/2 and i ^ k. Note 
that dU coincides with the Ising-DW profile (I2]),(l3b| for 
i = 1 and with the Bloch-DW profile (Il]),(|3al) for i = 2. 
In order to establish the magnetization dynamics in the 
limit of large negative values of time, I invert the motion 
of TO*^*-* utilizing the property m(*)(x, 0) — — m^'^*(— cc + 
2xgj,0), where Xq^ = xoi — ln(i/)/(2fci), and find 

1 _ „2j)i+ln(,y) i,i+ln(i/)/2 

-a ^ . . + i^,2 ^ (10) 



1 + e2r)i+ln(iy) 



1 -I- e2j7,+ln(i/) 



for w + ln(i^)/2 <C ??fc + ln(z^)/2 and i ^ k. Compar- 
ing ([8]) and ([TOl) . one sees that the result of the colhsion 
is a change of the sign of the real part of m^*^. Thus, one 
can think of the colliding DWs that they pass through 
each other with constant velocity and change their char- 
acter from the head-to-head into the tail-to-tail structure 
and vice versa. In fact, the result of the DW collision 
is their elastic reflection accompanied by change of the 
Bloch DW into Neel DW and vice versa, similarly to 
the result of field-driven or spontaneous collision of mag- 
netic DWs of the ID Landau-Lifshitz equation (I). Let 
me emphasize that taking the t — >■ —00 limit of .gi//i at 
^1(2) + ln(i^)/2 w leads to m differing from pUj) by the 
sign of its imaginary part, which confirms the irrelevance 
of (dl) in this limit. 

Above, I analyzed an infinite ID medium. Finite 
length of the domains leads to additional consequences 
of the DW collision because of the natural tendency of 
finite systems toward energy minimization (the energy 
cannot be defined for infinite systems) . The preferred di- 
rection of the DW motion before the collision (the growth 
or decrease of the bubbles) corresponds to the decrease 
of the Zeeman energy while the collision-induced motion 



in the reversed direction finishes when the decrease of 
the interaction energy of the (repulsing) DWs equals the 
Zeeman-energy increase. Thus, I predict the appearance 
of localized bubbles due to an external-field application 
in the parameter range /3i > 3/32. Their diameter is de- 
termined by the strength of the driving field. 

Approaching the Bloch-Ising transition point with /3i 
results in a decrease of (/3i — >• 3/3^ ^ v ^ 0+). Then, 
the imaginary parts of m, rh decrease and, eventually, 
vanish at /3i = 3/32. Simultaneously, the centers of both 
the DWs move away from each other due to the shift 
of the parameters ([7]) of (|6ap and (j6bp . [v ^ v + Ai/, 
kj -J> kj + Akj, Ak2 = 0). The shifts of the DW centers 
are equal to Axqj — \n{iy)/\2kj\ — \n{i'+Ai^)/\2kj+2Akj\. 
Finally, at /3i = 3/32, (Aiy = -v, k^ + Akj = ^2/32/J) 
the centers of both the DWs diverge. A domain pattern 
is expected to appear due to finite density of the mutu- 
ally interacting bubbles. Further decrease of /3i (down 
to —(32) results in divergence of the widths of the (Ising) 
DWs and, thus, in the divergence of the length of the DW 
interaction. An adiabatic approaching with /3i — >■ — /3^ 
is impossible in practice due to the instant widening of 
the DW. Because of the rapid increase of the overlap 
of DWs, one anticipates the appearance of domains of 
an unstable phase in our bistable medium in the regime 
—(32 < (3i < 3/32 (above the Bloch-Ising transition and 
below the bifurcation point). Then, an energy excess 
is being removed from the system via the propagation of 
fronts connecting domains of a stable phase with domains 
of an unstable phase (e.g. fronts between ferromagnetic 
and paramagnetic domains) . They always propagate into 
unstable state and below I call them phase fronts (29I ]. 



III. PHASE-FRONT COLLISIONS 

We can gain insight into the pattern formation near the 
criticality by considering the collision of fronts propagat- 
ing into unstable states. Let us study front solutions to 
([!]) assuming Im(m) = or, alternatively, to a simplified 
version of ([T]), 



a- 



dm 
dt 



J 



Om- ijm^ + -fH ^0, (11) 



where — /3i + (32, (the real solutions to both the equa- 
tions are the same). Equation (jlip is relevant to the clos- 
est vicinity of the critical point 9 = 0, where the order 
parameter takes real values. On the other hand, Im(m) 
vanishes when crossing the Bloch-Ising transition point 
/3i = 3/32; thus, real solutions are expected to qualita- 
tively describe the phase-front motion for the parameter 
range —(32 < (3i < 3/32. Since m ^ with 6* — >• 0, 
I introduce a renormalized field m = mj \fQ in order to 
determine dynamical properties of the phase fronts in the 
vicinity of the critical point. Substituting m by g//, one 



arrives at a trilinear (Hirota) form of ([TT|) : 

f[{~aDt + JDl)g-f + ^H/Vef] 

+g [{^JDl + e)f-f- fiOg^] = 0. 
For H = 0, an exact single- front solution 



^k{x — xai)+ckt 



/=1 



xoi )+ckt 



here 



ck = — , 
2a' 



(12) 
(13) 

(14) 



and a double-front solution to ([T 
1 



9 



(e 



k{x — xoi)-\-ckt 



-k{x — xo2)+':kt\ 



y _ 2 + gfe(2:-2:oi)+cfct _j_ Q-k{x-xo2)+ckt {15) 

were found by Nozaki and Bekki [SO^I. I mention that 
these authors have considered a generalization of pT|) ex- 
changing a, J, /i into complex coefficients (and allowing 
complex m), however, the real parts of these coefficients 
were found to determine generic dynamical properties of 
the phase fronts whereas the imaginary parts are respon- 
sible mainly for additional periodic structure of the phase 
domains [30], (ID spiral waves [sil). The solution (fT5|) 
describes a collision of the phase fronts at the borders of 
oppositely-oriented domains of the stable phase. Upon 
the collision both the fronts create the so called Nozaki- 
Bekki hole, a localized object which is similar to the Ising 
DW. 

For H — 0, following the derivation in Appendix B, we 
find an approximate double-front solution to p2)) of the 
form 

g _ ^ f ^k{x-xoi)+ckt ^-k(x-xo2)+ckt 
_^0Qk{x — xoi)+ckt^—k{x—xo2)+ckt^ 
J — }_ l^k{x~xoi)+ckt _|_ i^-k{x-xo2)+' 



kt 



+96" 



k= \ 



ck — — . 

2a 



(16) 



It describes two phase fronts at the borders of similarly- 
oriented domains of the stable phase since, contrary to 
the Nozaki-Bekki hole, all terms of g are of the same 
sign. They counterpropagate into an unstable area. Ac- 
cording to since i — 00 g f /y/JI, these fronts 
annihilate upon the collision. Above I have neglected the 
external field because it introduces a big complexity to 
the solution (Appendix B). However, the asymmetry of 
the Ginzburg-Landau function due to H ^ plays a role 
in the collision-induced formation of a dissipative struc- 
ture since it determines the choice of the preferred stable 
state [121 , [though, the strength of the field is limited by 
the condition j\H\ < 2{Q jifl'^ j p}'"^ ensuring that two 
stable homogeneous solutions to (ITT|) are of the opposite 
signs] . 




FIG. 1. 2D pattern structures in a bistable medium; (a) a 
bubble structure, (b) a lamellar (labyrinth) pattern, (c) stable 
(white) domains propagating into (grey) unstable state. In (c) 
a periodic (spiral) structure of the domains is neglected. 



IV. DISCUSSION 



Direct observation of the front motion in ferromag- 
netic and ferroelectric media is difficult because relevant 
time scales are very narrow. Thus, I refer to an exper- 
iment on a (bistable) ferrocyanide-iodate-sulfite chemi- 
cal reaction performed in a continuous-flow stirred tank 
reactor [l3, [ll|. When the ferrocyanide concentration 
is decreased, the morphology of the 2D chemical struc- 
ture changes. At a relatively-low concentration, irreg- 
ular stationary patterns of DWs (lamellae) appear and 
change into bubble structures with a further concen- 
tration decrease. The lamellae formation is accompa- 
nied by the property of the DW reflection upon their 
externally-driven collision (which is enforced by an in- 
tense irradiation). On increasing the critical parameter 
(the ferrocyanide concentration), the lamellae dissapear 
while domains become inhomogeneous (spiral structure 
of the domains appears) and their fronts annihilate upon 
the collision. I identify the transition between structures 
analogous to visualized in Figs. 1(a) and 1(b) with the 
decrease of the critical parameter and the transition 
between structures in Figs. 1(b) and 1(c) with crossing 
the Bloch- Ising point /3i = 3/32. Thus, the predicted dy- 
namical behavior of DWs of a generic (Ginzburg-Landau) 
model of bistable systems coincides with that observed 
in a chemical reaction below the critical point. Further 
increase of the ferrocyanide concentration leads to the 
appearance of wave-front reflection and self-replication 
of phase spots. These effects are related to differences in 
the diffusion coefficients of the chemical species; there- 
fore, they cannot be predicted within the Ginzburg- 
Landau model. According to calculations using a two- 
species (Gray-Scott) model of a chemical reaction, the 
self-replication regime relates to the monostability of the 
system [331]. 

I conclude that the coexistence of subcritical effects 
of front reflection and annihilation (accompanied by the 
creation of dissipative structures) can be explained by 
the appearance of domains of the unstable phase above 
the Bloch-Ising transition point. 
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Appendix A: Derivation of DW solutions 

One manages to find an exact single-DW solution to 
([T]) for 7^ in the case Im(m) — 0. It corresponds to 
an asymmetric Ising DW. In particular, inserting 



gi = ai- a2e~'', /i = 1 + e~ , 

rj ~ k{x — xq) — ujt, (Al) 

into (IT^ . one establishes the coefficients oi, 02, fc, w to be 
comphcated functions of H. For <C (/3i+/52)"^^^/a'^^^j 
{H is much weaker than the coercivity field), the effects 
of the DW asymmetry ai 7^ 02 and the dependence of 
the DW width 1/fc on H are neglig ible while re- 
produces ^ and (j3bp . Unfortunately, an exact solution 
describing a Bloch DW in an external field is unknown. 
Thus, I search for approximate solutions with Im(m) 7^ 
assuming they satisfy the following requirements; 1) they 
tend to the exact (single DW or double DW) solutions 
of Ref. HI] with H ^ 0, 2) they satisfy ^ at the cen- 
ters of the DWs (therefore, I neglect an expected space 
asymmetry of Re(m), Im(TO) with respect to the DW 
center in the presence of if 7^ 0). Within this approach, 
let us consider approximate single-DW solutions to the 
bilinearized Ginzburg-Landau system (j5a|) of the form 



31 = a(l-e2") + i6e'', /i = 1 + e 



2rj 



T] ~ k(x — X{)) — ujt, (A2) 



a{/3i +l32-X)+jH + e''i6(/3i - /?2 - A 



-He^''i6(/3i - /32- X + Jk^ - au) 

+c4''[-a(/3i+/32-A)+7i/] =0, 
-a^H + X + e2''(2aV - + 2A - 8Jk^) 

+e^''(-aV + A) =0. (A3) 



For X = f3i + /32, a — y//3i + /32/y^, one finds the sec- 
ond equation of (jA3P to be satisfied, (all coefficients of 
the LHS expansion in e'' vanish), if 6 = 2^//3i — 3(32/ ^/JI 
and I A: I = y^/^/j (the Bloch Dw) or 6 = and 
|fc| = y/PT+l^/V^J (the Ising DW), whereas, for uj = 



/JI-fH/{yW+^a), 
^31) simplifies to 



the LHS of the first equation of 



7iJ(l + e''ife 



,2^2 - e3"i6- 



+ 6^") 



We see that it vanishes with — or at 77 = (at 
the DW center). Therefore, the DW described by (TMj) 
propagates with a constant velocity c = u/k oc H. 

Let me mention that I do not perform a systematic 
Hirota construction connected to inserting consecutive 
elements of the expansion of gi , /i in e''^ , e^^ and solv- 
ing a series of relevant equations one by one. If their 
solutions were not exact, it would accumulate an error 
in the consecutive steps. In the present approach, the 
form of the solution ()A2p is fixed while I simply verify 
the coefficients a, b, uj. Similarly, searching for the ap- 
proximate double- DW solution, I insert (pa|) into (|5al) . 
The resulting equations take the form 



(1 + e4''i+'"'^)[a(/3i +j32~X)+ jH] + e''''ihv^'^{Pi - P2 - X + Jkf + auji) + c'^'^'2v{-fH - 2aauji) 
+c^'^-'2iy{jH - 2aauj2) + e^'^^ihv^'^iPi - ^2- X + Jkf- awi) + (e""! + c'^'^-')i^^[-a{l3i + /32 - A) + -yH] 
^g2^i+4r,2 2^(^ij + 2aauji) + c^''^^+'^'^''2v{-iH + 2aauj2) - ci'+^''Hbi^^/^{{l - i^)[(3i - (32 - X + Jkj + AJkj 
+au}i] + {1 + v){AJkik2 + 2auj2)} + c^'^^+^'^^2[a{l - u^){(3i +P2-X + AJkf + AJkl) + (1 + iy^){8aJkik2 
+-fH)] + e^^i+^^^i^Tyi/ZKi _ j,)[/3i _ _ A + Jkj + AJkl " a^i] + (1 + v){AJkik2 - 2a^2)} 

_e'Ji+4')2i^j,3/2^^^ - ^2 - A + J/c? + awi) - e3''i+4''2i6t/i/2(^^ ~ (32 - X + Jkl - auji) = 0, 
(1 -I- e^^ij/^ ^ + e4''i+4''^)(-a2/i + A) + (c^^i + e'^'i^+^'^^)v{2a^ + 2X- SJkf - fi) + e^^^+^^^2 

X [-(1 - iy^)8Jkik2 + (1 + ly^){-a'^^l + A - AJkl - '^J^D + vb'^^A + (e^*" + e^-'^^+^'^^)v{2a^ + 2X- SJk^) = 0(.A4) 



The second equation of (|A4p is fulfilled for A, a, b de- 
termined above and for ki, k2, i' of When H — 0, 
the first equation of (jA4ll is satisfied for this sort of pa- 
rameters according to [2^. In the vicinity of the DW 
centers in a state of well separated DWs; the regimes 
r?i + ln(j/)/2 « 0, e''=+''^(^)/2 ~ or 7^2 ln(7^)/2 « 0, 
gr,i-Hin(i.)/2 ^ -(-j^p system (jXI)) reproduces there- 



fore, the centers of both the DWs counterpropagate with 
constant velocities. 
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Appendix B: Derivation of double- front solution 



Looking for the approximate double-front solution 
(fTH]) . I define an error function as the absolute value of 
the LHS of (in for H = 



i:vTix,t) = \f[{-aDt + JDl)g-f] 
+g [{-JDl + 9)f-f~fi9g'] |. 



(Bl) 



Inserting 



g = h (e"! + e 



V2 _u ^pe'^i+m'j 



V1+V2 



(B2) 



where rji = k{x — xoi) + ckt, 772 = -~k(x — XQ2) + ckt, into 
(|B1I) , we take 1^9 = 6* in order to ensure faster divergence of 
the front centers in the limit — O"*" than the divergence 
of the front widths in this limit. We arrive at 



Err(a:,t) = | [e^^i+^^^ei^ + (e^^i+S''^ + e 



3jji+2j72 



+ (e' 



x6(l - + (c^'Ji + c2''2)6(2 - ack/e 
-Jk^/0) + 6^"^+^"^ 6*26(8 - Gfib^ ~ 2ack/0) 
+(e''i+2'?2 + e2''i+"^)06(7 - 3fib^ - 3ack/e 
+Jk'^/e) + e''i+''^6(5 - 4ack/e + QJk^/e) 
+(6^1 + e''^)6'-i6(l - ack/e + Jk^/O)]. (B3) 



For /i = l/T^i, fc = ^/e/2J, ck = 30 /{2a), all the coor- 
dinate X dependent terms of (|B3|) vanish, leading to 

Err(x,0 = Err(t) - |e2'7i+2';2 5/2 _ 2e''i+''^| 

= |e^'^'^*~2fe(a;oi -2:02)^2 _ 2g2cfct-fc(i:oi-2;o2) I (B4) 

Before the collision, when the fronts are well separated, 
Eq. ()12p [equivalent to Err(a;, t) — 0] is approximately 
fulfilled in the vicinity of the front centers [the regimes 
?7i +ln(6l) « 0, e''=+i"(^) « or r]2+\n{e) « 0, e''i+'"('') « 
0] as well as in the area of unstable phase separating the 
fronts (e''i+i'^(^) w 0, e''^+i"(^' w 0). Then 



1 



-l-e'' 



^ 1/6* + e''i -I- e''2 + 6'c''i+''; 

r 1 c"i+'"<''' 



(B5) 



_ I yj7 l+o-il+i"!") 

^ S 1 c';2+ln(«) 

I 7^ l + 0'J2+ln(f) 



V2 



ln(6'),e''2+'"(^) - 0, 
ln(6l),e''i+i"(^) ~ 0. 



Assuming the motion to start at t = 0, and choosing 
the reference frame such that xq2 = —xqi, the front 
centers are initially at the points x = ±xi = ±(xoi — 
\n{9)y/2J/9) and they meet at the point a; = after 
the time t — x\/c. Let us notice that up to the time 
t = t + the function Err(t) is a decreasing one and 
Err(t ) = 0. We conclude that the approximate solution 
([T6l) satisfies well ([T2l) at the beginning of the motion (at 
t = 0) and the accuracy of approximating the solution 
by the functions (|16l) increases with time up to t = t 
which corresponds to the moment of finishing the front 
collision (the moment when the fronts of the width 1/fc 
finish their passing through each other with the relative 
velocity 2c). 

The application of a similar analysis to the field in- 
duced collision of DWs below the Bloch-Ising transition, 
studying ((5a|) separately from (jSbl) . does not work. Uti- 
lizing the argumentation (of Sec. II) based on the time- 
reversal symmetry is an alternative. 
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